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Abstract 

A  propeller  lifting-surface  design  and  analysis  program  is  improved  upon  by  implementing 
enhancements  in  the  source  distribution  calculation  to  represent  the  blade  thickness.  It  is 
recognized  that  the  present  method  of  setting  the  source  line  distribution  representing  blade 
thickness  (currently  based  on  linearized  slender-body  theory  for  an  isolated  foil  section)  may 
introduce  significant  errors.  This  is  the  case  for  propulsors  with  a  combination  of  a  large 
thickness/ chord  ratio  (blockage  effect)  and  numerous  blades  (cascade  effect). 

A  source  panel  (area)  method  was  developed  to  more  accurately  model  these  effects.  This 
method  uses  the  lattice  structure  of  the  current  PBD-14  code  from  which  to  compute  the  source- 
induced  velocity  factors  between  the  blades,  hub,  and  duct,  if  present.  Using  the  method 
of  images  allows  the  hub  and  duct  to  be  modeled  as  panel  images  from  the  blade  panels. 
The  source-induced  effects  of  the  whole  propulsor  are  accounted  for  by  using  a  panel  method 
to  obtain  a  source  distribution  along  the  mean  camber  surface  of  the  blade.  Invoking  the 
kinematic  boundary  condition  on  the  true  blade  suction  and  pressure  surfaces  solves  this  system 
of  linear  equations,  which  represent  the  blade  thickness  distribution.  This  robust  formulation 
assigns  source  strengths  more  accurately  over  a  much  larger  range  of  thickness/chord  ratios 
and  increasing  numbers  of  blades,  as  evidenced  by  a  more  accurate  velocity  streamline  trace 
representation  of  the  actual  pressure  and  suction  side  surfaces  of  the  blade.  Experimental 
validation  is  demonstrated  for  open  and  ducted  fiow  stators. 

Thesis  Supervisor:  Justin  E.  Kerwin 
Title:  Professor  of  Naval  Architecture 
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Nomenclature 


Mathematical  Notation 


A 

area  of  quadrilateral  panel  (flat  2:  —  plane) 

BB 

blade-to-blade  influence  functions 

DB 

duct-to-blade  influence  functions 

dR 

radial  difference  between  source  element  endpoints 

ds 

differential  component  of  local  coordinate  position 

dS 

differential  component  on  the  blade  surface 

dx 

axial  distance  between  control  points  in  a  chordwise  strip 

dt 

6t 

differential  thickness  between  control  points  in  a  chordwise  strip 

ductrat 

^bldelem 

^ductelem 

duct  element  imaged  source  strength  adjustment  factor 

D 

propeller  diameter 

hi,h2 

ds 

metric  coefficient 

HB 

hub-to-blade  influence  functions 

HIF 

horseshoe  influence  function  matrix 

hubrat 

L>bld^lem 

^huhelem 

hub  element  imaged  source  strength  adjustment  factor 

1 

pressure  side  of  the  blade  (lower  surface) 

Rbldelem 

length  of  blade  source  element 

Lduclelem 

length  of  blade  duct  element 

Rhubelem 

length  of  blade  hub  element 

total  length  of  a  source  element 

m 

mean  camber  surface  of  the  blade 

n 

normal  vector  on  a  surface 

P 

surface  coordinate 

P 

arbitrary  point  in  space 

R 

radial  coordinate 

R 

X(S1,S2)  -x(s'i,5'2) 

differential  position  vector 

Sl,S2 

local  coordinates 

S 

blades  surface 

SIP 

source  influence  function  matrix 
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t 


u 

Ui, 

Us 


Ve// 

^ind 


^indfiuc 
^  total 

Vis 

Vas 

w 

X 

X 

y 

z 

+/- 

a 

$ 

7 

r 

n 

O-A 

a 

(^L 

^^src 

^^snk 

YjAsrc 


wake  surface  created  by  horseshoe 
tangential  vector  on  a  surface 
H  velocity  associated  with  the  x-coordinate 
inflow  velocity  to  the  propeller 
component  of  inflow  velocity 
||  velocity  associated  with  the  y-coordinate 
effective  velocity 
total  induced  velocity 
induced  circumferential  mean  velocity 
fluctuating  component  of  induced  velocity 
total  velocity 

component  of  velocity  due  to  vortices 
component  of  velocity  due  to  sources 
^  velocity  associated  with  the  z-coordinate 
surface  coordinate  vector  for  x,  y,  z 
axial  coordinate 
vertical  coordinate 
z-coordinate 

upper/lower  blade  surface  (suction/pressure 
angle  of  attack  of  propeller  blade 
potential 

vortex  element  strength  (per  unit  length) 

blade  circulation 

blade  rotation  rate  (^) 

source  panel  distribution  (per  unit  area) 

source  strength  distribution  (2D  sense) 

source  line  distribution  (per  unit  length) 

sum  of  positive  source  strengths 

sum  of  negative  source  strengths  (sinks) 

sum  of  panel  areas  due  to  source  strengths 

sum  of  panel  areas  due  to  sink  strengths 
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Abbreviations 


CMV 

horseshoe 

BANS 


circumferential  mean  velocity 
a  vortex  structure  composed  of;  a  spanwise  vortex, 
its  shed  chordwise  vortices  on  the  blade  and  in  the 
transition  wake,  and  the  ultimate  wake  helices. 
Reynolds- Averaged  Navier  Stokes 
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Chapter  1 


Introduction 


1.1  Overview 

The  goal  of  this  research  was  to  create  a  more  accurate  model  which  describes  the  true  lifting 
surface  of  a  propeller  blade  system.  This  system  includes  the  interactions  between  the  propeller 
blades  and  the  hub,  as  well  as  the  duct,  if  present.  Since  many  current  propeller  system  design 
trends  include  the  use  of  a  multiple-staged  ducted  propulsors,  accurately  modelling  the  effects  of 
induced  velocities  between  the  different  components  of  the  propulsor  are  paramount  to  perform 
a  good  analysis.  There  are  several  numerical  tools  and  techniques  which  accomplish  this  task 
sufficiently  [1]: 

1.  Lifting-Line  Methods 

2.  Vortex-Lattice  Lifting-Surface  Methods 

3.  Boundary  Integral  Equation  Methods 

4.  Axisymmetric  and  Three  Dimensional  Reynolds  Averaged  Navier  Stokes  (RANS)  Flow 
Solvers. 

The  vortex-lattice  lifting-surface  technique  provides  an  accurate  and  robust  method  to  examine 
the  design  [2]  and  analysis  [3]  of  complex  marine  propulsors  efficiently.  Kerwin  et  al.  [2] 
[3]  developed  such  a  lifting-surface  technique  which  has  been  improved  upon  [4]  [5]  since  its 
introduction. 
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The  vortex-lattice  lifting-surface  technique  is  a  potential  flow  method  which  does  not  fully 
describe  what  occurs  within  the  propulsor  system,  like  full  viscous  flow  models.  However, 
Kerwin  et  al.  [4]  developed  a  coupling  procedure  which  uses  the  potential  flow  lifting-surface 
technique  in  conjunction  with  a  RANS  viscous  flow  solver.  This  coupled  technique  allows  for 
multiple  blade  row  analysis  while  providing  an  accurate  model  of  the  physical  behavior  of  the 
system. 

The  linearized  slender-body  theory  representation  does  not  accurately  model  the  flow  over 
the  propeller  blades  because  it  under  predicts  the  effects  of  the  thickness  distribution  of  the 
propeller  blades.  It  also  does  not  account  for  the  blade- to-blade,  blade-to-hub,  and  blade-to- 
duct  interactions.  Therefore,  as  more  blades  and  stages  are  added  to  the  propeller  system, 
the  linearized  slender-body  theory  method  is  unable  to  accurately  predict  the  cascade  and 
blockage  effects  created.  This  deficiency  in  the  current  vortex-lattice  method  is  the  motivation 
to  implement  an  enhanced  thickness  model  in  an  effort  to  better  approximate  the  flow  patterns 
through  the  system. 

1.2  Background 

1.2.1  Vortex-Lattice  Lifting-Surface  Analysis  Tools 

A  comprehensive  explanation  of  the  vortex-lattice  lifting-surface  theory  and  methods  for  the 
design  and  analysis  of  marine  propulsors  is  documented  extensively  [2],  [3],  [6],  and  [7].  Propeller 
Blade  Design  (hereafter  designated  PBD-I4)  developed  by  Kerwin^  et  al.  was  used  for  this 
research.  In  PBD-14^  a  lattice  geometry  of  discrete  vortex  segments  is  placed  along  the  blade 
mean  camber  surface.  Control  points  are  placed  within  each  grid  lattice.  Blade  thickness  is 
represented  using  a  linearize  slender-body  theory  to  obtain  a  source  strength  distribution  at 
the  lattice  segments  in  a  strip- wise  sense.  Each  lattice  segment  has  known  source  strength 
and  the  vorticity  strength  is  solved  for,  such  that  the  kinematic  boundary  condition  is  satisfied 
along  the  mean  camber  surface,  namely  the  total  velocity  must  be  tangent  to  the  surface  at 
each  control  point.  A  trailing  wake  system  is  grown  from  the  trailing  edge  of  the  blades.  See 
Figure  1-1  for  the  vortex-lattice  discretization  of  the  propeller  geometry. 

^Justin  E.  Kerwin,  Profesor  of  Naval  Architecture,  Ocean  Engineering  Department,  MIT. 
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Figure  1-1:  Vortex  lattice  discretization  of  a  propeller  geometry  (4119)  where  the  symbols 
represent  control  points. 

1.2.2  Propeller  Blade  Thickness  Modeling 

It  is  recognized  that  the  present  method  of  setting  the  strength  of  the  source  distribution 
representing  the  blade  thickness  (based  on  linearized  slender-body  theory  for  an  isolated  foil) 
may  introduce  significant  errors  in  the  case  of  propulsors  with  a  combination  of  a  large  thick¬ 
ness/chord  ratio  and  numerous  blades.  The  questions  to  be  answered  are: 

1.  Is  the  thickness-induced  velocity  correct  in  PBD-14^  There  are  three  parts  to  this  prob¬ 
lem: 


•  isolated  two-dimensional  foil, 

•  two-dimensional  cascade  effects  (emulations  of  multi-blade  flow)  and 

•  three-dimensional  effects. 


15 


2.  How  does  one  accurately  model  these  effects  in  a  marine  propulsor  system?  Two  methods 
will  be  discussed  below: 

•  linearized  slender-body  theory  and 

•  panel  method  for  a  source  distribution  along  the  mean  camber  surface. 

Due  to  efforts  to  more  accurately  assign  source  strengths  to  the  vortex-lattice  structure 
(i.e.  correctly  model  cascade  and  blockage  effects  in  propulsors),  enhancements  were  sought  in 
the  current  methodology  of  the  source  distribution  to  represent  the  blade  thickness  distribu¬ 
tion.  The  current  version  of  PBD-I4  uses  linearized  slender-body  theory  to  predict  the  source 
strength  necessary  to  satisfy  a  pseudo  kinematic  boundary  condition  at  the  designated  control 
points.  Regardless  of  the  blade  geometry,  the  kinematic  boundary  condition  is  invoked  for  flow 
over  a  thickened  foil,  assuming  zero  camber  and  loading.  It  performs  a  systematic  assignment 
of  source  strengths,  one  chord  wise  strip  at  a  time  using  equations  1.1  and  1.2: 

/X2  rx2 

adx  =  Vtotai  /  -j-dx  (1.1) 

X\  J  Xi 

which  simplifies  to 

O'L—^totalSt  (1.2) 

where 


(j 

o-L 


^  total 

6t 


dx 


the  source  strength  distribution 

the  source  strength  distribution  per  unit  length 

the  total  inflow  velocity  in  the  presence  of  an  operating  propeller 

the  increase  in  blade  thickness  between  lattice  control  points 

along  a  given  chordwise  strip 

the  axial  distance  between  control  points  being  compared  to 
determine  the  source  strength 
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X 

Figure  1-2:  Comparison  of  desired  thickness  form  and  source-induced  thickness  form  of  a  Bi¬ 
convex  blade  section  obtained  by  streamline  tracing  for  current  PBD-I4  code. 

Therefore,  the  strength  of  each  source  line  is  determined  using  the  relationship  between  the 
total  inflow  velocity  and  the  differential  thickness  of  the  blade  at  some  position  (x)  downstream 
from  the  leading  edge  of  the  blade  section.  Figure  1-2  shows  an  extreme  test  case  for  an  infinite 
pitch  stator  in  uniform  flow  with  zero  circulation  (F)  and  a  thickness/chord  ratio  of  27  percent. 
The  desired  thickness  form  (solid  line)  can  be  compared  with  the  actual  thickness  form  (vector 
symbols)  obtained  by  streamline  tracing.  The  source-induced  thickness  form  of  the  section  is 
about  20  percent  too  thin,  thus  indicating  the  source  strength  is  too  small. 

Source  Line  Method 

In  this  method,  the  source  distribution  is  determined  for  the  vortex-lattice  structure  by  modeling 
the  complete  blade  system  (all  blades,  the  hub,  and  the  duct,  if  present).  This  is  a  vortex-lattice 
method  which  distributes  the  source  line  elements  on  the  mean  camber  surface  of  a  propeller 
blade.  Then,  the  source  distribution  for  the  entire  blade  system  is  determined  by  satisfying 
the  kinematic  boundary  condition  at  control  points  which  are  on  the  actual  blade  surface 
(pressure  and  suction  sides).  The  result  is  a  source  strength  per  unit  length  at  the  vortex/source 
elements.  This  yields  extremely  accurate  results  for  the  foil  shape  shown  in  Figure  1-2,  but  gets 
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progressively  inaccurate  once  the  thickness/chord  ratio  falls  below  ten  percent. 

Source  Panel  Method 

A  alternative  and  more  robust  method,  determines  the  source  distribution  for  the  complete 
blade  grid  structure  by  fully  modeling  the  entire  blade  system  (all  blades,  the  hub,  and  the 
duct,  if  present)  with  source  panels.  This  is  a  panel  method  which  distributes  the  source  panels 
on  the  mean  camber  surface  of  a  propeller  blade.  Then,  the  source  distribution  for  the  entire 
blade  system  is  determined  by  satisfying  the  kinematic  boundary  condition  at  control  points 
which  are  on  the  actual  blade  surface  (pressure  and  suction  sides).  The  result  is  a  source 
strength  per  unit  area  which  can  be  integrated  over  the  blade  grid  structure  to  obtain  a  similar 
lattice  source  strength  per  unit  length  as  above.  This  yields  extremely  accurate  results  for  the 
foil  shape  shown  in  Figure  1~2  for  all  ranges  of  the  thickness/chord  ratio. 

1.2.3  Hub  £ind  Duct  Modeling 

The  hub  is  represented  by  a  similar  vortex-lattice  geometry  which  is  appropriately  placed  using 
the  method  of  images.  The  duct,  if  present,  is  developed  in  the  same  manner.  The  constant- 
strength  source  line  elements  in  the  hub  and  duct  are  ratioed  from  the  blade  source  element 
strengths  by  comparing  the  lengths  of  the  blade  elements  with  their  associated  imaged  elements 
in  the  hub  and  duct.  Figure  1-3  shows  a  hub  and  duct  representation  in  PBD-14- 

1.3  Objective 

This  thesis  provides  a  description  of  the  advanced  thickness  modeling  techniques  implemented 
into  the  lifting-surface  propeller  blade  design  and  analysis  code  PBD14-3  [8].  This  enhancement 
extends  the  capability  of  current  lifting-surface  theory  by  properly  accounting  for  thickness 
effects  (cascading  and  blockage)  in  this  model  structure  more  completely.  An  overview  of  the 
process  is  presented  first,  followed  by  more  detailed  descriptions  of  the  thickness  modeling 
process.  Finally,  validation  examples  are  presented. 
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Figure  1-3:  Blade,  hub  and  duct  representation  in  PBD-I4  for  Ka-455  propulsor. 
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Chapter  2 


Lifting-Surface  Design  and  Analysis 
Theory 


2.1  Overview 

In  the  past,  lifting-surface  theory  has  provided  accurate  results  by  using  slender-body  theory 
thereby  approximating  higher-order  thickness  and  higher-order  viscous  effects  for  conventional 
propulsors  as  equivalent.  It  has  been  shown  that  in  conventional  thickness  and  loading  distri¬ 
butions,  the  decrease  in  lift  due  to  viscous  effects  is  offset  by  the  increase  in  lift  due  to  blade 
thickness  [2] .  Due  to  the  advances  in  blade  section  design,  the  higher-order  viscous  and  higher- 
order  thickness  effects  do  not  necessarily  cancel.  Hence,  these  factors  must  be  accounted  for  in 
the  design  and  analysis  of  marine  propulsors. 

Lifting-surface  propeller  design  and  analysis  places  a  vortex-lattice  structure  on  the  mean 
camber  surface  of  the  blade.  The  presence  of  thickness  is  accounted  for  by  the  distribution  of 
sources  and  sink  lines  collocated  with  the  vortex  elements.  There  is  a  local  force  as  a  result 
of  the  velocities  of  a  source,  influencing  its  neighboring  vortex  elements.  An  inter-blade  effect 
occurs  due  to  the  blockage  of  the  modeled  thickness  of  the  blade.  There  is  also  a  Lagally  force 
that  arises  from  the  force  on  the  sources  themselves.  Panel  method  codes  model  the  thickness 
of  the  blade,  so  these  effects  are  addressed.  However,  current  lifting-surface  codes  do  not  have 
this  ability  [3]  [9].  Also,  blade  thickness  results  in  a  blockage  to  the  axisymmetric  flow  which 
affects  the  flow  distribution  past  the  propulsor  and  duct,  if  present.  This  factor  must  also  be 
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taken  into  consideration  if  propulsor  performance  is  to  be  predicted  more  accurately  [10] . 

2.2  Vortex-Lattice  Method 

The  vortex-lattice  method  of  solving  the  lifting-surface  propeller  blade  problem  can  be  catego¬ 
rized  into  four  separate  steps: 

1.  A  continuous  vortex/source  sheet,  located  on  the  blade’s  mean  camber  surface,  is  linearly 
represented  by  a  discretized  geometry  for  the  propeller  blade. 

2.  Determine  the  influence  coefficients  for  the  vortex/source  distributions. 

3.  Apply  the  kinematic  boundary  condition  which  stipulates  that  the  flow  on  the  blade 
surface  is  only  in  the  tangential  direction. 

4.  Solve  the  system  of  linear  equations  developed  from  the  steps  outlined  above. 

2.2, 1  Geometry  Discretization 

Numerically,  circulation  and  thickness  distributions  are  discretized  using  a  lattice  structure 
of  constant-strength,  straight-line  vortex/source  elements  that  are  aligned  in  the  chordwise 
and  spanwise  directions.  These  vortex/source  elements  are  collocated  along  the  blade’s  mean 
camber  surface.  PBD-I4  uses  a  vortex-lattice  technique  which  incorporates  a  uniform  B-spline 
representation  of  the  blade’s  mean  camber  surface  [2].  This  representation  of  the  blade  allows 
for  efficient  interrogation  of  the  B-spline  control  polygon  to  discretize  the  blade’s  surface  [6]. 
Furthermore,  various  lattice  structure  spacing  options  as  available  (i.e.  cosine  and/uniform 
spacing)  as  described  by  Greeley  and  Kerwin  [3].  Kelvin’s  Theorem  states  that  vortex  lines 
cannot  end  anywhere  in  the  flow  field,  therefore  trailing  vortices  convect  downstream.  The 
bound  element  with  its  associated  convected  trailers  form  a  horseshoe  vortex  which  extends 
beyond  the  trailing  edge  of  the  propeller  blade  [3] .  The  trailing  vortices  must  be  aligned  with 
the  flow  in  order  to  be  force-free. 
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2.2.2  Influence  Coefficients 


Lifting-surface  methodology  solves  the  propulsor  problem  by  first  developing  an  influence  ma¬ 
trix,  This  matrix  uses  the  induced  velocity  produced  by  each  vortex/source- lattice 

segment  as  the  influence  (assumed  unit  strength  of  one)  upon  each  of  the  blade  control  points. 
Therefore,  the  influence  coefficients  depend  solely  on  the  geometry.  When  this  influence  matrix 
is  multiplied  by  the  actual  vortex  segment  strength,  the  induced  velocity  at  every  control  point  is 
known.  Section  2.2.4  more  fully  explains  the  influence  coeflBcients  purpose  in  the  lifting-surface 
methodology. 

2*2.3  Kinematic  Boundary  Condition 

In  order  to  correctly  describe  the  flow  through  the  propeller  blade  system,  one  must  implement 
the  appropriate  boundary  condition  on  the  rigid  body.  The  kinematic  boundary  condition 
stipulates  that  no  flow  will  penetrate  the  surface  of  the  blade  (i.e.  the  total  velocity  flow  is 
tangential  at  the  surface).  Another  way  of  stating  the  kinematic  boundary  condition  is  that  the 
total  velocity  normal  to  the  grid  surface  must  be  zero.  This  condition  is  specified  at  a  number 
of  control  points  on  the  mean  camber  surface  which  is  equal  to  the  number  of  vortex/source 
elements  representing  the  blade  lattice  structure  and  is  shown  in  equation  2.1  below; 

^totai  •  n  ==  0.  (2.1) 

Section  2.2.4  more  fully  explains  how  the  kinematic  boundary  condition  is  utilized  to  solve  the 
lifting-surface  problem. 

2.2.4  Solution  Procedure 

The  lifting  surface  can  be  defined  by  a  curvilinear  system  of  coordinates  (si,S2)  on  the  blade 
surface  as  shown  in  Figure  2-1.  The  surface  coordinates  are  defined  as 

X(S1,S2)  =  W5l,S2),y(si,S2),-2(Sl,S2)]  (2.2) 

and  the  metric  coefficients  are  given  as 
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(2.3) 


ds 


and 


h2  = 


dx 

dS2 


(2.4) 


Figure  2-1:  Vortex  lattice  discretization  of  a  propeller  geometry  where  the  symbols  represent 
control  points  and  si  and  S2  are  orthogonal  vectors  on  the  blade  surface. 


The  kinematic  boundary  condition  is  given  as 


[7(gi,4)  X  R]  -n 


hih2ds[ds'2 


-[C^m(5l,  S2)  +  +  ^ind(si,  Sa)]  •  n 


(2.5) 


where 
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R  =  x(si,S2) -x(si,4) 

n  =  vector  normal  to  S  at  point  (51,52)- 

Equation  2.5  specifies  that  there  is  no  fiow  through  the  blade  surface  S',  considering  a 
continuous  sheet  of  vorticity  (7)  on  the  blade  surface  and  in  the  wake  (Sw)  in  the  presence  of 
an  incoming  flow  (Uin)  and  rotational  velocity  {QR).  Since  the  thickness  distribution  for  the 
blade  is  known,  the  source-induced  velocities  (Vindtotai )  ^^7  computed  and  included  on  the 
right  hand  side  of  the  equation  [10]. 

A  lifting-surface  code  is  used  for  either  a  design  of  a  new  propulsor  or  an  analysis  of  an 
existing  propulsor.  In  each  case,  equation  2.5  applies,  but  the  known  versus  unknowns  vary. 
For  propulsor  design,  a  desired  radial  and  chordwise  loading  distribution  is  prescribed.  Then, 
the  blade  shape  is  manipulated  until  equation  2.5  is  satisfied.  For  propulsor  analysis,  the 
blade  shape  is  prescribed.  Then,  the  resulting  circulation  is  solved.  This  process  summation  is 
simplified  to  matrix  equation  2.6: 


[HIFj-[r,]  =  -[V,.n]  (2.6) 

where 

HIFij  =  matrix  of  horseshoe  influence  coefficients.  The  subscripts  denote  the  influence  of 
the  horseshoe  vortex  on  the  control  point.  Each  element  represents  the 
velocity  component  normal  to  S',  induced  at  a  control  point  due  to  a  horseshoe 
vortex  of  unit  strength. 

Fj  =  blade  circulation  scalar  composed  of  individual  vortex  strengths  of  the 
vortex. 

V j  •  fi  =  velocity  component  normal  to  the  blade  mean  camber  surface  at  the  control 
point. 
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Solving  this  matrix  system  of  linear  equations  yields  equation  2.7: 


[r,]  =  -[V,-n][HIFj-i  (2.7) 

Blade  thickness  effects  are  added  by  placing  discrete  source  lines  coincident  with  the  chord- 
wise  blade  vortex  elements.  The  resulting  propulsor  forces  due  to  the  lifting  surface  and  thick¬ 
ness  are  calculated  from  the  Kutta-Joukowski  and  Lagally  theorems,  respectively.  A  leading 
edge  suction  force  and  Lighthill  pressure  distribution  correction  are  applied  to  those  forces,  and 
the  propulsor’s  sectional  viscous  drag  is  calculated  based  on  either  stripwise  two-dimensional 
empirical  drag  coefficients  or  a  stripwise  two-dimensional  integral  boundary  layer  calculation. 

2.3  Propulsor  Inflow 

The  accuracy  of  modeling  the  propulsor  is  dependent  upon  correctly  representing  the  propeller 
inflow.  The  three  concepts  for  propeller  inflow  are: 

Definition  1  Nominal  inflow:  Velocity  field  at  the  plane  of  the  propeller  when  the  propeller  is 
not  operating. 

Definition  2  Total  inflow:  Velocity  field  at  the  plane  of  the  propeller  when  the  propeller  is 
operating. 

Definition  3  Effective  inflow:  Total  inflow  without  the  effects  of  the  propeller  induced  veloci¬ 
ties. 

In  propeller  design  and  analysis,  the  lifting-surface  methodology  deals  with  the  flow  resulting 
from  the  operation  of  the  propeller  in  a  specified  axisymmetric  effective  inflow.  Therefore,  the 
total  velocity  relation  at  the  blade  surface  in  the  axial,  radial  and  tangential  directions  can  be 
decomposed  into  the  following  three  components 


^ total  —  ^eff  “b  ^indcmv  “b  ^indfiuc 


(2.8) 


where 


25 


^ total 


^eff 


LioLLcrnv 
tTld  flue 


=  the  total  velocity  in  the  presence  of  an  operating  propeller  at  a  given  point  on 
the  blade. 

=  the  effective  wake  behind  the  ship. 

=  the  circumferential  mean  velocity  induced  due  to  propeller  blades  and  wakes. 
=  the  fluctuating  component  of  induced  velocity  due  to  propeller  blades. 


Therefore,  the  total  velocity  is  a  combination  of  the  effective  inflow  and  the  net  propeller 
blade  induced  velocities.  In  the  propeller  blade  problem,  the  total  velocity  is  used  to  solve  the 
kinematic  boundary  condition  as  shown  in  equation  2.1. 
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Chapter  3 


Numerical  Implementation  of  the 
Enhanced  Blade  Thickness  Model 


3.1  Overview 

In  the  current  version  of  PBD-14^  propeller  lifting-surface  design  and  analysis  program  code, 
the  thickness  distribution  is  developed  in  the  following  manner  using  linearized  slender-body 
theory: 

1.  For  a  chordwise  strip,  determine  the  thickness  distribution  from  the  B-splined  input  file 
parameters  and  the  change  in  thickness  (St)  between  the  vortex-lattice  structure  control 
points  (dx). 

2.  Use  the  flow  field’s  circumferential  mean  velocity  (CMV)  and  rotational  velocity  compo¬ 
nent  {Q.R)  to  determine  the  blade’s  local  velocity  at  the  control  points  (along  the  mean 
camber  surface). 

3.  Solve  for  the  source/sink  line  elements  constant-strength  {crLij)  necessary  to  create  the 
correct  6t  between  the  control  points  on  that  spanwise  strip.  Equation  3.1  shows  that  the 
source  strength  per  unit  length  is  a  function  of  inflow  velocity,  the  change  in  the  thickness 
distribution  and  the  differential  length  between  control  points: 

~  (2*^) 
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4.  Adjust  each  source  line  element  strength  for  “slant.”  This  effect  increases  the  strength  of 
the  source  line  element  when  it  is  not  perpendicularly  aligned  with  the  inflow  by  a  factor 
of  .  Here,  Ltotai  is  the  total  length  of  the  element  and  dR  is  the  radial  difference 
between  element  endpoints. 

5.  Correct  the  sink  strengths  to  ensure  each  chordwise  strip  of  source/sink  line  elements 
sums  to  zero.  This  effect  increases  each  sink  element  strength  by  a  factor  of  thus 
ensuring  that  the  <Tsrc  +  E  ^snk  =  0. 

This  methodology  was  flawed  for  cases  where  the  blade  design  had  numerous  blades  and/or 
a  high  thickness/chord  ratio  (>  10%).  This  scheme  did  not  account  for  blade-on-blade,  hub- 
on-blade,  or  duct-on-blade  effects  in  the  source  line  strength  determination.  Therefore,  as  the 
source  lines  of  each  modeled  propeller  section  and  its  associated  components  get  closer  together, 
these  effects  increase  in  importance  (i.e.  far  field  effects  begin  to  move  into  the  near  field  range 
to  increase  the  influence  of  these  elements  on  one  another) . 

Because  of  the  shortcomings  of  this  methodology,  more  accurate  methods  of  accounting  for 
blade,  hub  and  duct  interactions  to  more  fully  represent  the  propulsor  were  investigated.  Two 
promising  techniques  were  developed,  implemented  and  tested  in  PBD-I4 : 

1.  Source  Line  Method 

2.  Source  Panel  Method. 

The  source  panel  method  is  an  improvement  upon  the  work  completed  in  the  source  line 
method.  This  work  was  undertaken  to  develop  a  more  robust  representation  of  all  the  effects 
that  determine  the  blade  thickness  distribution.  For  example,  the  source  line  method  accounted 
for  the  cascade  effect  (increasing  the  number  of  blades  creates  larger  blade-to- blade  influences). 
It  also  was  able  to  properly  represent  the  blade  thickness  distribution  for  large  thickness/chord 
ratios.  Both  of  these  efforts  were  improvements  over  the  shortcomings  of  the  previous  PBD-I4 
methodology  of  assigning  source  strengths.  However,  the  source  line  method  over  predicted  the 
blade  thickness  distribution  when  the  thickness/chord  ratio  decreased  to  less  than  ten  percent. 
This  effect  got  worse  as  the  thickness/chord  ratio  decreased.  The  source  panel  method  corrects 
the  shortcomings  of  the  source  line  method,  so  that  it  accurately  represents  the  blade  thickness 
distribution  over  the  broadest  range  of  thickness/chord  ratios  and  numbers  of  blades. 
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3.2  Source  Line  Method 


This  thesis  seeks  to  find  improvements  in  PBD-I4  to  more  accurately  represent  the  blade 
thickness  distribution  over  the  blade  lattice  structure  in  the  presence  of  numerous  blades  or 
large  thickness/chord  ratios  led  to  the  development  of  the  source  line  methodology.  An  outline 
of  this  method  follows: 

1.  Discretize  the  blade  geometry  on  the  mean  camber  surface  so  that  the  vortex  and  source 
elements  are  collocated. 

2.  Develop  each  chordwise  blade  thickness  distribution  from  the  PBD-14  input  file  ~  values 
shown  in  Appendix  B. 

3.  Fully  develop  the  blade’s  actual  outer  surface  for  each  spanwise  position  at  each  vortex 
element  endpoint  that  lies  on  the  mean  camber  surface.  This  is  accomplished  by  creating 
the  pressure  and  suction  side  surfaces  from  the  actual  thickness  distribution  of  the  blade 
section  at  the  mean  the  camber  surface  vortex/source  elements.  The  mean  camber  surface 
lattice  endpoints  were  moved  in  the  normal  direction  by  the  appropriate  6t  for  the  specified 
point.  Figures  3-1  and  3-2  show  the  discretized  geometry  of  a  propeller  blade  with  the 
pressure  and  suction  side  surfaces  included. 

4.  Develop  a  complete  set  of  control  points  and  a  normal  to  the  control  points  on  the  pressure 
and  suction  side  of  the  actual  blade  surface.  This  is  done  in  the  same  fashion  as  the  mean 
camber  surface  calculations  in  previous  versions  of  PBD-14- 

5.  Develop  the  hub  and  duct  imaged  lattice  structure.  The  hub  and  duct  images  lie  along 

the  mean  camber  surface  of  the  inner  and  outermost  spanwise  set  of  endpoints,  thus 
having  the  same  pitch  of  the  blade  at  the  hub  and  duct  intersection  points.  Also,  the 
hub  and  duct  imaged  source  strengths  are  adjusted  by  a  factor  of  and 

respectively. 

6.  Determine  the  total  inflow  velocity  at  the  control  points  on  the  blade’s  mean  camber 
surface. 
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7.  Develop  a  system  of  linear  equations  in  matrix  form,  similar  to  equation  2.6,  to  solve 
for  the  source  strengths  associated  with  each  of  the  source-induced  function  coefficients 
[SIFij] .  This  is  done  while  satisfying  the  kinematic  boundary  condition  iy total  •  n  =  0) 
at  the  control  points  on  the  actual  pressure  and  suction  side  surfaces  of  the  blade.  The  set 
of  matrix  equations  is  increased  to  add  a  further  condition  that  the  Gsrc  +  Yl,  <^snk  =  0 
for  each  chordwise  strip  on  the  blaide. 

This  set  of  matrix  equations  is  shown  in  a  simplified  notation  in  equation  3.2: 

ISIFjjl  ■  K1  =  -|Vj  ■  nl  (3.2) 

where 

SlFij  “  matrix  of  source  influence  coefficients.  The  subscripts  denote  the  influence  of 
the  source  element  on  the  control  point.  Each  element  represents  the 
velocity  component  normal  to  5,  induced  at  a  control  point  due  to  a  source 
element  of  unit  strength. 

=  constant  strength  source  composed  of  individual  source  strengths  of  the 
source  elemnts. 

Vj  •  ri  =  velocity  component  normal  to  the  blade  mean  camber  surface  at  the 
control  point. 

Solving  this  matrix  system  of  linear  equations  yields  equation  3.3 

K]  =  -[Vra][SIF,J-^  (3.3) 

In  order  to  determine  what  truly  affects  the  blade  thickness  distribution  problem,  a  decom¬ 
position  of  source  and  vortex  velocities  at  the  pressure,  suction  and  mean  camber  surfaces  was 
necessary.  Figure  3-3  shows  a  simple  chordwise  blade  representation  for  the  velocity  effects  on 
the  blade  thickness  distribution.  The  proof  follows. 
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Figure  3-1:  Xy-coordinates  for  blade  thickness  discretization  on  the  pressure  and  suction  sides. 
Definitions  for  the  nomenclature  in  Figure  3-3  is: 


-f  = 
m  = 

n  = 
t  = 

(S1,S2)  = 


{Us^,U,,)  = 


(Vas^Vas^)  = 


pressure  surface  (lower)  of  the  blade, 
suction  surface  (upper)  of  the  blade, 
mean  camber  surface  (middle)  of  the  blade, 
normal  for  pressure  and  suction  side  surfaces, 
tangent  for  pressure  and  suction  side  surfaces, 
local  coordinates  in  normal  (1)  and  tangential  (2) 
directions  on  the  mean  camber  surface, 
component  of  inflow  velocity  in  the  normal  (1) 
and  tangential  (2)  directions, 
component  of  velocity  due  to  sources, 
component  of  velocity  due  to  vortices. 
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Figure  3-2:  XZ-coordinates  for  blade  thickness  discretization  on  the  pressure  and  suction  sides. 

In  order  to  meet  the  kinematic  boundary  condition  on  both  the  pressure  and  suction  side 
surfaces  of  the  blade,  equations  3.4  and  3.5  must  hold  true. 


(3.4) 

(3.5) 

where 

n+  =  —  n_  and  =  t_  exactly, 
and  for  small  thickness  distributions, 

Vlsi+  ^ 

From  substitution  and  adding  equations  3.4  and  3.5,  the  following  result  is  obtained 


{Usx  +V(^si  +^1^751)  ^  +  (^«2  +  V^S2  +  ^7^2)  * 


+ 


=  0 


{Usi  +Vas^  +  ^7si)  II  +  {Us2  +  Vo's2  +  ^7*2)  ^ 


=  0 


{Vas,^-Vas,.) 


n+  + 


Us2  + 


{V(Ts2^+V(Ts2-)  ,  (^7.2+  +^7.2-) 


+ 


t_  =0. 


(3.6) 
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Figure  3-3:  Chordwise  strip  of  a  blade  section  to  represent  the  components  which  affect  blade 
thickness  distribution. 


Assuming  a  small  thickness  distribution,  Figure  3-4  shows  a  simplified  mean  camber  line  rep¬ 
resentation  of  T^7s2- 

Since  F752,u  =  +  2  ~  ^^S2,m  ~  2>  averaging  the  effects  of  the  vortex 

induced  velocities  from  the  pressure  and  suction  surfaces  of  the  blade  at  the  mean  camber 

surface  produces  ^7^2  =  2  ^  ~  ^^S2,m'  H^^ice,  equation  3.6  can  be  simplified  further 


n+  -f- 


Us2  + 


(F, 


<^S2+  +  F(Ts2_) 


+  F7. 


S2 


I'-l- 


(3.7) 


Furthermore,  for  a  symmetric,  two  dimensional  foil,  ^7^^  =  0-  Also,  if  the  foil  has  a  moderate 
camber,  like  most  propeller  blade  sections,  ^7^^  fa  0  and  can  be  considered  as  a  part  of  the 
CMV  inflow.  This  assumes  that  is  a  small  percentage  of  total  CMV  inflow  which  is  made 
up  of  many  components.  Therefore,  the  component  of  velocity  can  be  moved  to  the  right 
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Figure  3-4:  Mean  camber  line  representation  of  vortex-induced  velocities  averaging  to 

hand  side  of  equation  2.6  and  iterated  upon  to  get  a  final  solution. 

Finally,  and  Vas^_  on  the  pressure  and  suction  side  surfaces  of  the  blade  can  be 

expressed  as  a  dot  product  of  the  source-induced  influence  coefficients  [SIF^  and  the  unknown 
source  line  strength  [(Tl^Y  Thus,  equation  3.7  simplifies  to  the  matrix  equation  3.2. 

To  check  this  reasoning,  a  symmetrical,  two-dimensional  foil  is  placed  at  an  arbitrary  angle 
of  attack  (a)  as  shown  in  Figure  3-5.  Since  =  C/i^iSina,  Us^_  —  C/i^sina  and  ^42+  — 
Us2-  ”  iJmCOsa,  then  substituting  these  parameters  into  equation  3.7,  yields  the  expected 
result  that  the  blade  thickness  distribution  is  only  affected  by  Uin  cos  a  or 


Vasj^n-^  +  {Uin  cos  a  +  Vas2)t-^  =  0,  (3.8) 

Once  source-induced  coefficients  are  developed  for  the  pressure  and  suction  side  surfaces 
of  the  blade,  they  are  averaged  to  get  a  representative  [SIF^  j]  coefficient  matrix  for  the  mean 
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Figure  3-5:  Two-dimensional  blade  section  used  to  check  what  attributes  affect  blade  thickness 
distribution. 

camber  surface.  The  mathematical  formulation  is  solved  using  a  least  squares  methodology 
which  already  existed  in  the  PBD-I4  code.  In  this  solution  formulation,  the  hub-to-blade  and 
duct-to-blade  imaged  influence  coefficients  are  added  directly  to  the  blade-to-blade  influence 
coefficients.  This  is  explained  more  fully  in  Section  3.4.2.  This  is  a  more  robust  and  accurate 
method  than  the  current  PBD-I4  code  since  the  kinematic  boundary  condition  is  satisfled  on 
the  actual  blade  surface  at  the  new  set  of  control  points  developed  on  this  surface. 

3.3  Source  Panel  Method 

A  second  method  to  accurately  model  the  source  distributions  over  the  blade’s  mean  camber 
surface  was  developed  after  the  shortcomings  of  the  source  line  method  were  found.  While 
the  source  line  method  accurately  represents  propulsors  with  numerous  blades  and  large  thick- 
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ness/chord  ratios,  it  fails  to  correctly  represent  thinner  blades  where  the  thickness/chord  ratio 
approaches  zero.  In  fact,  it  begins  to  over  predict  the  source  line  strengths  and  hence  the  blade 
thickness  distribution  when  the  thickness/chord  ratio  becomes  less  than  about  ten  percent. 

Because  source  line  elements  are  a  set  of  discrete  concentrations  of  constant-strength  sources, 
they  tend  to  average  the  effects  between  control  point  positions,  where  the  kinematic  boundary 
condition  is  invoked  for  equation  3.9.  This  becomes  more  difficult  for  the  system  of  linear 
equations  as  the  blade  thickness  distribution  approaches  zero  (i.e.  matrix  equation  3.2  becomes 
ill-conditioned  and  obtaining  a  solution  by  the  least  squares  method  is  impossible).  Therefore, 
the  next  logical  step  to  obtain  a  less  discrete  concentration  of  the  constant  strength  source 
distributions,  but  still  accurately  represent  the  blade  thickness  distribution  is  to  move  from  line 
sources  {(Tl)  to  area  sources  (a a)  or  panels.  This  source  panel  methodology  seeks  to  obtain 
results  with  the  same  accuracy  as  the  source  line  method  at  large  thickness/chord  ratios  and 
numerous  blades.  However,  it  also  seeks  to  improve  upon  the  linearized  slender-body  theory 
approximation  results  of  the  current  PBD-I4  code  as  thickness/chord  ratio  approaches  zero. 
This  source  panel  method  is  outlined  below: 

1.  Discretize  the  entire  B-spline  blade  geometry  on  the  mean  camber  surface,  to  include  both 
the  leading  and  trailing  edges.  This  means  that  one  extra  control  point  must  be  created 
at  the  leading  edge. 

2.  Develop  each  chordwise  blade  thickness  distribution  from  the  PBD-14  input  file  ^  values 
shown  in  Appendix  B. 

3.  Fully  develop  the  blade’s  actual  outer  surface  for  each  spanwise  position  at  each  vortex 
element  endpoint  that  lies  on  the  mean  camber  surface.  This  is  accomplished  by  creating 
the  pressure  and  suction  side  surfaces  from  the  actual  thickness  distribution  of  the  blade 
section  at  the  mean  the  camber  surface  vortex/source  elements.  The  mean  camber  surface 
lattice  endpoints  were  moved  in  the  normal  direction  by  the  appropriate  6t  for  the  specified 
point.  Figures  3-1  and  3-2  show  the  discretized  geometry  of  a  propeller  blade  with  the 
pressure  and  suction  side  surfaces  included. 

4.  Develop  a  complete  set  of  control  points  and  a  normal  to  the  control  points  on  the  pressure 
and  suction  side  of  the  actual  blade  surface.  This  is  done  in  the  same  fashion  as  the  mean 
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camber  surface  calculations  in  previous  versions  of  PBD-I4.  However,  the  leading  and 
trailing  edge  control  points  are  moved  off  the  actual  edge  of  the  mean  camber  surface  to 
a  position  between  the  lattice  structure  elements. 

5.  Develop  the  hub  and  duct  imaged  panel  structure.  The  hub  and  duct  images  lie  along  the 

mean  camber  surface  of  the  inner  and  outermost  spanwise  set  of  endpoints,  thus  having 
the  same  pitch  of  the  blade  at  the  hub  and  duct  intersection  points.  However,  the  hub  and 
duct  imaged  source  strengths  no  longer  require  an  adjustment  by  a  factor  of  and 

respectively.  This  is  due  to  the  robust  nature  of  the  panel  methodology  imposed 
with  this  technique. 

6.  Determine  the  total  inflow  velocity  at  the  control  points  on  the  blade^s  mean  camber 
surface. 

7.  Develop  a  system  of  linear  equations  in  matrix  form,  similar  to  equation  3.2,  to  solve 
for  the  source  strengths  associated  with  each  of  the  source-induced  function  coefficients 
[SIF^^j] .  This  is  done  while  satisfying  the  kinematic  boundary  condition  (Vtotai  *  n  =  0) 
at  the  control  points  on  the  actual  pressure  and  suction  side  surfaces  of  the  blade.  The  set 
of  matrix  equations  is  increased  to  add  a  further  condition  that  the  X)  ^src  +  ^snk  ~  0 
for  each  chord  wise  strip. 

This  set  of  matrix  equations  is  shown  in  a  simplified  notation  in  equation  3.9: 


[SIFj.[a^J  =  -[V,--n].  (3.9) 

Similar  to  the  reasoning  discussed  in  Section  3.2,  the  same  conclusions  can  be  drawn  about 
the  velocity  contributions  to  the  blade  thickness  distribution.  The  following  section  will  cover 
the  modifications  required  to  the  PBD-I4  code  to  implement  this  methodology. 

New  influence  function  subroutines  have  to  be  coded  with  the  equations  found  in  Appendix 
A.  In  the  three-dimensional  case,  the  discretization  has  two  parts: 

•  discretize  the  geometry 

•  determine  the  singularity  element  distribution. 
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If  these  elements  are  represented  by  a  B-spline  net,  both  geometry  and  singularity  strength, 
then  the  first  order  approximation  to  the  surface  can  be  defined  as  a  quadrilateral  panel^ 
with  a  constant-strength  source  singularity.  Since  a  vortex-lattice  mean  camber  surface  is 
represented  as  a  lattice  structure  of  rectilinear  panels,^  a  conversion  from  a  rectilinear  surface 
to  a  quadrilateral  surface  is  performed  [11].  From  Appendix  A,  the  potential  at  an  arbitrary 
point  P(x,  y,  z)  due  to  this  quadrilateral  element  is 


Hx,y,z)  =  — 


47r  Js  + 


(3.10) 


and  the  velocity  components  can  be  obtained  by  differentiating  the  velocity  potential: 


(3.11) 


This  diflFerentiation  from  equation  3.11  results  in  source-induced  velocity  coefficients  used  to 
obtain  blade- to-blade,  hub-to-blade,  and  duct-to-blade  influence  functions. 


3.4  Integration  of  the  Source  Panel  Model  in  PBD-I4 

3.4.1  Blade  Lattice  Modification 

The  current  methodology  uses  a  vortex-lattice  structure  which  does  not  include  the  actual 
leading  edge  section  (B-spline  leading  edge  endpoints  to  the  first  vortex/source  element).  The 
new  methodology  is  to  use  the  entire  cubic  B-spline  blade  structure,  leading  edge  to  trailing 
edge,  to  describe  the  source  distribution  lines.  This  allowed  for  a  symmetric  distribution  of 
source  panels  for  cosine  spacing  in  the  chordwise  direction  along  the  mean  camber  surface. 
However,  the  number  of  source  panels  increases  by  one  more  than  vortex  elements.  Control 
points  were  placed  as  before  with  two  exceptions: 

1.  a  leading  edge  control  point  was  added  and  placed  in  a  similar  fashion  as  the  other  control 
points  (about  the  middle  of  the  lattice  rectilinear  panel);  and 

quadrilateral  panel  is  a  flat  surface  with  four  straight  sides. 

rectilinear  panel  has  straight  but  not  necessarily  flat  sides  that  can  be  twisted. 
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2.  the  trailing  edge  control  point  was  moved  and  placed  in  a  similar  fashion  of  the  other 
control  points,  vice  right  on  the  trailing  edge. 

These  two  exceptions  were  necessary  so  that  the  control  points  did  not  lie  on  the  edge  of  a 
source  panel  and  create  a  singularity  for  the  system  of  linear  equations,  thus  forming  an  ill- 
conditioned  matrix.  Similarly,  the  hub  and  duct  images  were  developed  from  source  panels  vice 
source  line  elements. 

To  further  enhance  the  model  and  truly  satisfy  the  kinematic  boundary  condition,  the  lattice 
points  for  the  entire  mean  camber  surface  were  moved  to  both  the  pressure  and  suction  side 
surfaces  of  the  blade.  From  these  new  surfaces,  control  points  were  added  in  the  same  fashion 
as  the  current  PBD-I4  code  did  on  the  mean  camber  surface.  Figures  3-1  and  3-2  show  the 
fully  developed  blade  lattice  structure. 

Lastly,  the  rectilinear  panels  which  form  the  entire  mean  camber  surface  set  of  source  panels 
had  to  be  transformed  into  quadrilateral  panels  which  approximate  the  same  constant-strength 
source  distribution.  This  transformation  was  completed  so  that  the  code  in  reference  [11], 
which  solves  for  aA,  could  be  utilized.  This  transformation  was  accomplished  using  a  set  of 
subroutines  for  PSF10.3  code.  These  subroutines  input  the  rectilinear  panel  ends  consisting 
of  four  endpoint  coordinates.  Then,  the  subroutines  return  the  quadrilateral  panel  endpoints 
which  lie  on  a  constant  z  —  plane  and  a  set  of  coordinates  for  the  actual  center  of  the  original 
rectilinear  panel.  Therefore,  the  orientation  between  the  local  and  global  coordinate  system  is 
maintained. 

3.4.2  Adjustments  to  the  Solution  Procedure 

The  solution  procedure  for  constant-strength  source  area  distribution  is  completely  different 
from  the  current  PBD-I4  code.  It  now  accounts  for  blade-to-blade,  hub-to-blade  and  duet¬ 
to-blade  influence  interactions  which  affects  induced  velocities  produced  within  the  propulsor. 
These  new  influence  coefficients  are  the  averaged  values  of  the  pressure  and  suction  side  [SIF^ 
matrix  which  represents  the  mean  camber  surface  [SIF^^j]  matrix.  The  solution  is  solved  as  a 
system  of  linear  equations  via  a  least  squares  method  vice  a  strict,  two-dimensional,  linearized, 
slender-body  theory.  By  incorporating  the  effects  of  the  other  blades,  the  hub  and  the  duct, 
if  present,  the  model  more  accurately  represents  cascade  and  blockage  effects  which  affect  the 
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blade  thickness  distribution.  Therefore,  the  propulsor  blade  representation  is  more  accurate  in 
a  larger  variety  of  propulsor  types. 

Since  panels  are  now  utilized  to  obtain  the  constant-strength  source  distributions,  there 
is  a  discrepancy  in  how  PBD-I4.  must  account  for  this  source  distribution.  One  method  is 
to  integrate  the  source  panel  strengths  in  the  chordwise  direction.  This  integration  technique 
would  lump  the  effects  of  the  panels  into  source  line  elements  which  are  collocated  with  the 
vortex  elements.  While  this  method  is  still  an  improvement  over  the  source  line  method  and 
linearized  slender-body  theory,  it  still  diminishes  the  robust  nature  of  the  method. 

Therefore,  a  separate  velocity  matrix  was  developed  from  the  new  known  source  strength 
distribution  for  the  blade,  hub  and  duct  interactions  at  the  old  control  point  positions  on  the 
mean  camber  surface. 


[Vindj [a Aj].  (3.12) 

Equation  3.12  shows  the  velocity  matrix  that  will  be  imported  into  the  [HIFij]  system  of 
linear  equations,  as  a  known  velocity  vector  of  source-induced  effects  on  the  right-hand  side  of 
the  equation.  This  velocity  vector  will  be  calculated  using  the  source  panel  influences  on  the 
specified  control  points  vice  the  source  line  element  influences.  Then,  this  modified  version  of 
the  [HIFjj]  system  of  equations  will  be  solved  to  obtain  [Fj]. 

3.4.3  Adjustments  to  Hub  suid  Duct  Modeling 

The  method  of  images  is  used  in  PBD-I4  to  represent  the  hub  and  duct.  This  provides  a  means 
to  implement  boundary  conditions  necessary  in  potential  flow  theory.  For  singularities  on  rigid 
boundaries,  potential  flow  introduces  another  singularity  into  the  flow  field  which  mirrors  the 
original  singularity.  The  vortex  element  images  model  the  hub  and  duct  accurately  and  are 
satisfactory  from  the  standpoint  of  meeting  the  kinematic  boundary  condition.  However,  the 
source  element  images  are  less  accurate.  The  resulting  imaged-lattice  structure  is  exhibited  in 
Figure  1-3.  The  imaged  source  panels  for  the  hub  and  duct  are  created  in  the  same  manner  as 
the  blade  source  panels. 
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Role  of  the  Hub  and  Duct 


Adding  the  hub  and  duct  source  influence  functions  to  the  [SIFij]  matrix  system  expands 
the  complexity  of  appropriately  accounting  the  influences  on  each  control  point,  but  does  not 
complicate  the  actual  system  of  linear  equations  that  must  be  solved.  Equation  3.13  presents 
the  expanded  influence  function  coefficient  for  a  generalized  control  point  position: 


Source  Panels 

Control  I  SIFBBiJ+SIFHBi.j  +  SIFDBij 
Points  ^ 


■[(taj]  =  -[Vj  ■  n]. 


(3.13) 


The  above  equation  appears  complex,  however,  it  is  merely  an  extension  of  the  blade-only 
formulation  (see  equation  3.9). 

Hub  and  duct  surfaces  at  the  blade  endpoints  are  streamlines  of  the  flow,  where  the  compo¬ 
nent  of  y^totai  normal  to  the  blade  surface  is  zero.  The  hub  and  duct  source  panel  influences  on 
the  blade  vortex-lattice  system  induces  zero  CMV  normal  to  the  blade  surface.  This  has  been 
exhibited  by  a  field-point  velocity  calculation  routine  at  the  control  points  on  the  blade  surface 
[6],  as  well  as  by  the  velocity  streamline  traces  developed  in  the  next  chapter. 
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Chapter  4 


Validation 


In  this  chapter,  two  methods  of  evaluating  stators  are  compared.  First,  an  infinite-pitch,  con¬ 
stant  ^  stator  with  a  hub  and  duct  in  a  uniform  inflow  is  evaluated.  The  original  PBD-I4 
scheme  will  be  compared  to  the  source  panel  method  for  a  variety  of  thickness  distributions  at 
different  angles  of  attack  to  exhibit  the  differences  in  how  each  method  represents  of  source- 
induced  effects.  Second,  an  infinite-pitch,  constant  ^  stator  in  a  uniform  inflow  is  evaluated 
at  the  hub  and  duct  to  ensure  their  effects  are  accounted  for  properly.  Finally,  the  number  of 
blades  will  be  increased  to  eleven  and  the  effects  analyzed  at  the  hub.  Since  the  most  signifi¬ 
cant  source-induced  factors  occur  at  this  point,  then  cascade  effect  analysis  will  be  performed 
here.  The  source  panel  method  will  exhibit  its  robustness  when  it  uses  varying  blade  thickness 
distributions  for  two  different  blade  section  geometries. 


4.1  The  Infinite-Pitch,  Constant  ^  Stator 

This  test  case  shows  that  the  source  panel  method  is  more  robust  than  the  linearized  slender- 
body  theory  of  the  original  PBD-I4  code  for  a  larger  range  of  thickness/chord  ratios  and  number 
of  blades.  Results  vary  with  vortex-lattice  grid  density  and  therefore  all  comparisons  will  be 
made  with  a  15  x  15  grid.  Table  4.1  shows  how  the  parameters  for  the  stator  will  be  varied  in 
this  test  case. 
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Blade  Section  Shape 

Angle  of  Attack  (a) 

t/D 

thickness  /  chord 

NACA66  modified 

0° 

0.020 

26.67% 

Biconvex 

1° 

0.015 

20.00% 

5° 

0.010 

13.33% 

0.002 

2.67% 

Table  4.1:  Infinite-pitch,  constant  t/D  stator  parameters 


X 

y 

z 

R 

0.00 

0.20 

0.00 

0.20 

0.05 

0.50 

0.00 

0.50 

0.10 

0.80 

0.00 

0.80 

0.15 

1.00 

0.00 

1.00 

Table  4.2:  Infinite-pitch,  constant  t/D  stator  with  a  zero  degree  angle  of  attack. 


4,1. 1  Initial  Comparison  for  Large  ^  Values  at  a  Zero  Angle  of  Attack 

Comparison  of  current  PBD-14  and  source  panel  method  for  the  blade  section  shapes  listed 
in  Table  4.1  is  exhibited  in  Figures  4-1  and  4-2.  This  set  of  figures  shows  how  the  linearized 
slender-body  theory  approximations  for  blade  thickness  distribution  are  inadequate  when  the 
blade  section  becomes  thicker.  Figures  4-3  and  4-4  show  the  results  for  the  source  panel  method 
which  has  less  than  one  percent  error  associated  with  the  blade  thickness  distribution.  Because 
of  the  straightforward  shape  (no  camber,  as  seen  in  Table  4.2)  and  uniform  inflow,  this  set  of 
stators  does  not  produce  any  circulation  or  forces  which  was  verified  by  the  PBDOUT.SGR  file. 
The  blade  lattice  structure  is  a  15  X  15  grid  with  hub  and  duct  images  (7)  implemented.  The 
dividing  streamline  is  shown  so  that  the  modeled  thickness  distribution  (velocity  streamlines) 
can  be  compared  with  the  actual  thickness  distribution  on  the  pressure  and  suction  side  surfaces 
of  the  blade.  The  cut  is  made  midspan  on  the  blade. 


4.1.2  Comparisons  at  a  One  Degree  Angle  of  Attack 

The  next  step  was  to  place  the  blades  at  an  angle  of  attack  greater  than  zero.  This  would 
test  equation  3.8  results  that  predicted  the  tangential  velocity  component  on  the  mean  camber 
surface  {Uin  cos  a)  is  the  only  affect  the  blade  thickness  distribution.  Figures  4-5  through  4-6 
exhibit  the  continued  trend  for  the  original  PBD-14  code  to  under  predict  the  blade  thickness 
distribution  at  the  larger  thickness/chord  ratios.  Also,  the  results  of  equation  3.8  are  verified  in 
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Figure  4-1:  Linearized  slender-body  theory  velocity  streamline  results  for  a  NACA66  modified 
blade  section  with  ^  =  0.020  and  a  =  0°  in  a  uniform  inflow. 


X 

y 

z 

R 

0.00000000 

0.20000000 

0.00000000 

0.20 

0.04999239 

0.19999810 

0.00087262 

0.50 

0.09998477 

0.19999239 

0.00174524 

0.80 

0.14997715 

0.19998287 

0.00261786 

1.00 

Table  4.3:  Infinite-pitch,  constant  t/D  stator  with  an  one  degree  angle  of  attack. 

Figures  4-7  and  4-8.  Because  of  the  angle  of  attack  (Table  4.3  shows  the  stator’s  B-spline  net 
design  for  this  a)  and  uniform  inflow,  this  set  of  stators  produces  some  circulation  and  forces 
which  was  verified  by  the  PBDOUT.SGR  file. 

4.1.3  Comparisons  at  a  Five  Degree  Angle  of  Attack 

The  next  step  was  to  place  the  blades  at  an  angle  of  attack  greater  than  one.  This  would  test 
equation  3.8  results  further.  Figures  4-9  and  4-10  exhibit  the  continued  trend  for  the  original 
PBD-U  code  to  under  predict  the  blade  thickness  distribution  at  the  larger  thickness/chord 
ratios.  Also,  the  results  of  equation  3.8  are  well  supported  in  Figures  4-11  and  4-12,  even  though 
the  streamlines  do  not  exactly  match  the  actual  thickness  distribution  shape.  However,  this 
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Figure  4-2:  Linearized  slender-body  theory  velocity  streamline  results  for  a  Biconvex  blade 
section  with  =  0.020  and  a  =  0°  in  a  uniform  inflow. 


X 

y 

z 

R 

0.00000000 

0.20000000 

0.00000000 

0.20 

0.04980973 

0.19985742 

0.00435779 

0.50 

0.09961947 

0.19971046 

0.00871557 

0.80 

0.14942920 

0.19957226 

0.01307336 

1.00 

Table  4.4:  Infinite-pitch,  constant  t/D  stator  with  a  five  degree  angle  of  attack. 


series  of  figures  definitely  shows  the  improvements  that  the  source  panel  method  makes  to  the 
overall  blade  thickness  distribution  representation  in  many  different  scenarios.  Because  of  the 
angle  of  attack  (Table  4.4  shows  the  stator’s  B-spline  net  design  for  this  a)  and  uniform  inflow, 
this  set  of  stators  produces  some  circulation  and  forces  which  was  verified  by  the  PBDOUT.SGR 
file. 


4.1.4  Varying  the  Thickness/ Chord  Ratio  at  a  Zero  Angle  of  Attack 

To  show  the  robustness  of  the  source  panel  method  over  the  Unearized  slender-body  theory, 
different  thickness/chord  ratios  are  compared  in  Figures  4-13  through  4-16,  including  the  Figures 
4-2  and  4-4  (thickness/chord  ratio  of  26.67%).  As  as  the  thickness/chord  ratio  approaches  zero. 
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Figure  4-3:  Source  panel  method  velocity  streamline  results  for  a  NACA66  modified  blade 
section  with  ^  =  0.020  and  a  =  0°  in  a  uniform  inflow. 

the  source  panel  method  and  linearized  slender-body  theory  begin  to  merge  to  the  same  blade 
thickness  distribution,  as  expected.  Note  that  as  the  blade  thickness  distribution  approaches 
zero,  it  is  better  to  increase  the  grid  density  to  obtain  a  finer  source  and  vortex  strength 
distribution  over  the  blade.  In  this  manner,  the  blade  representation  is  more  accurate,  but 
computational  time  increases.  For  some  cases  it  may  be  necessary  to  increase  the  grid  density 
to  25  X  25  or  even  35  x  35. 

4.1.5  Hub  cOid  Duct  Interfaces  and  the  Cascade  Effect 

Finally,  the  next  example  is  just  a  check  that  the  source  panel  method  with  hub  and  duct 
modeling  yields  the  correct  blade  thickness  distribution  at  these  intersecting  points  on  the 
model.  Figures  4-17  and  4-18  exhibit  the  robustness  of  the  source  panel  method  to  include 
all  effects  (blade-to-blade,  hub-to-blade,  and  duct-to-blade)  of  the  sourced-induced  functions 
accurately.  A  special  case  of  the  hub  interface  with  the  blades  is  shown  in  Figure  4-19  where 
the  number  of  blades  increases  to  eleven.  The  source  distribution  increases  to  appropriately 
represent  the  cascade  effect  in  the  model.  Note  that  the  source  panel  method  still  exhibits 
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Figure  4-4:  Source  panel  method  velocity  streamline  results  for  a  Biconvex  blade  section  with 
=  0.020  and  a  =  0°  in  a  uniform  inflow. 

its  robust  nature  to  represent  the  blade  thickness  distribution  at  the  hub  and  blade  interface 
accurately. 
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Figure  4-5:  Linearized  slender-body  theory  velocity  streamline  results  for  a  NACA66  modified 
blade  section  with  ^  =  0.020  and  a  =  1°  in  a  uniform  inflow. 


Figure  4-6:  Linearized  slender-body  theory  velocity  streamline  results  for  a  NACA66  modified 
blade  section  with  ^  =  0.020  and  a  =  1°  in  a  uniform  inflow. 
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Figure  4-7:  Source  panel  method  velocity  streamline  results  for  a  NACA66  modified  blade 
section  with  ^  =  0.020  and  a  =  1°  in  a  uniform  inflow. 


Figure  4-8:  Source  panel  method  velocity  streamline  results  for  a  Biconvex  blade  section  with 
=  0.020  and  a  =  1°  in  a  uniform  inflow. 
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Figure  4-9:  Linearized  slender-body  theory  velocity  streamline  results  for  a  NACA66  modified 
blade  section  with  ^  =  0.020  and  a  =  5°  in  a  uniform  inflow. 


Figure  4-10:  Linearized  slender-body  theory  velocity  streamline  results  for  a  Biconvex  blade 
section  with  ^  =  0.020  and  a  =  5°  in  a  uniform  inflow. 
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Figure  4-11:  Source  panel  method  velocity  streamline  results  for  a  NACA66  modified  blade 
section  with  =  0.020  and  a  5°  in  a  uniform  infiow. 


Figure  4-12:  Source  panel  method  velocity  streamline  results  for  a  Biconvex  blade  section  with 
^  =  0.020  and  a  =  5^  in  a  uniform  infiow. 
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Figure  4-13:  Linearized  slender-body  theory  velocity  streamline  results  for  a  Biconvex  blade 
section  with  ^  =  0.010  and  o;  =  0°  in  a  uniform  inflow. 


Figure  4-14:  Source  panel  method  velocity  streamline  results  for  a  Biconvex  blade  section  with 
^  =  0.010  and  a  =  0°  in  a  uniform  inflow. 
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Figure  4-15:  Linearized  slender-body  theory  velocity  streamline  results  for  a  Biconvex  blade 
section  with  ^  =  0.002  and  a  —  0°  in  a  uniform  inflow. 
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Figure  4-17:  Source  panel  method  velocity  streamline  results  for  a  Biconvex  blade  section  with 
=  0.015  and  a  =  0°  in  a  uniform  inflow  at  the  hub/blade  interface. 


Figure  4-18:  Source  panel  method  velocity  streamline  results  for  a  Biconvex  blade  section  with 
^  =  0.015  and  a  =  0°  in  a  uniform  inflow  at  the  duct/blade  interface. 
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Figure  4-19:  Source  panel  method  velocity  streamline  results  for  a  Biconvex  blade  section  with 
=  0.015  and  a  =  0°  in  a  uniform  inflow  at  the  hub/blade  interface.  This  streamtrace  is  for 
the  key  blade  of  an  11-blade  stator  system. 
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Chapter  5 


Conclusion 


5.1  Source  Panel  Method  Results 

In  this  thesis,  several  enhancements  to  the  blade  thickness  representation  for  lifting-surface 
design  and  analysis  for  marine  propulsors  were  discussed.  Optimization  over  the  whole  propulsor 
allowed  for: 

•  more  accurate  representation  of  the  blade  thickness  distribution.  This  was  accomplished 
by  incorporating  the  effects  of  blade-to~ blade  interaction  (blockage  effect),  as  well  as  the 
hub-to-blade  and  duct-to-blade  influence  effects. 

•  still  able  to  use  the  original  PBD-I4  code  structure  to  obtain  solutions  since  the  blade 
thickness  distribution  is  represented  by  constant-strength  source  panels  whose  effect  is  in¬ 
corporated  into  an  overall  source-induced  velocity  term  on  the  right-hand  side  of  equation 
2.6. 

•  more  accurate  results  for  larger  ranges  of  thickness/chord  ratios  and  numbers  of  blades 
(cascade  effect). 

The  weakness  of  the  source  panel  method  is  that  the  computational  time  increases  when  a 
solution  of  linear  equations  is  sought  to  determine  the  source  strength  distribution.  Also,  as  a 
finer  grid  density  is  used  to  obtain  more  accurate  results,  the  computational  time  increases. 

However,  this  panel  method  enabled  the  kinematic  boundary  condition  to  be  satisfied  on 
the  actual  pressure  and  suction  side  surfaces  of  the  blade.  This  ensures  that  Y total  on  the  blade 
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surface  is  tangential.  The  hub  and  duct  source  influences,  as  well  as  other  blades  influences, 
were  shown  in  the  validation  results  of  Chapter  4  for  the  following  conditions: 

•  varying  thickness/chord  ratios, 

•  vary  angles  of  attack  and 

•  two  distinct  blade  section  geometries. 

The  velocity  streamline  traces  showed  that  the  source  panel  method  is  more  accurate  and 
robust  over  a  the  specified  range  of  thickness/chord  ratios,  angles  of  attack  and  blade  section 
geometries. 

5.2  Recommendations  for  Future  Enhancements 

This  thesis  implemented  a  panel  method  into  the  PBD-I4  code  along  the  mean  camber  surface 
of  the  vortex-lattice  structure.  Some  recommendations  for  further  improvements  and  future 
work  in  this  area  include: 

•  integrate  the  source  panel  method  with  the  vortex  lattice  hub  and  duct  formulation  devel¬ 
oped  in  PBD-14-4  coding.  This  would  allow  for  a  more  accurate  representation  of  the  hub 
and  duct,  as  well  as  represent  their  associated  affects  on  the  blade  thickness  distribution. 

•  develop  a  more  accurate  field  point  velocity  grid  calculation  for  the  entire  propulsor  vice 
just  a  chordwise  strip  to  capture  all  three-dimensional  effects  visually. 

•  carry  the  source-induced  velocity  vector  from  the  panel  method  solution  to  calculate  a 
more  accurate  force  and  pressure  distribution  over  the  actual  pressure  and  suction  side 
surfaces  of  the  blade. 
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Appendix  A 


Constant-Strength  Quadrilateral 
Source  Distribution  Formulation 
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TO 


Y 


x-j; 


(xi,yi,  0) 
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A.l  Near  Field  Formulation 

Execution  of  the  integration  within  the  area  bounded  by  the  four  straight  lines  requires  a  lengthy 
process.  The  details  of  integration  are  shown  in  reference  [12]  and  the  results  obtained  for  the 
potential  of  a  planar  element  becomes 


1  ri-{-r2+di2 

ri+r2~di2 


,  r2+r3+d23 

’"2+^3— C^23 

.  In  ^3+^4+^34 

^"3 +’"4—^34 


^x,y,z)  = 


where 


tan-i 
+  tan~^ 


.  In  r4+ri+d4i 
r4+ri— ^41 


+  tan-^ 


+  tan" 


1  /  7712363-/13^ 

V  ^^3  J 

1  f  m4iei-^\ 


di2  =  J{X2~  xiY  +  {y2  -  y\Y 


^23  =  \/{xs-  X2)‘^  +  (ys  -  y2)^ 


d^A  =  J{xa-  x^Y  +  (^4  -  y3)2 


dAi  =  J{xi-  xaY  +  {yi  -  y^Y 


y2-yi 

X2  -  Xi 


ys  -  j/2 

Xs  ”  X2 
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(A.IO) 


m34  = 


X4  -xa 


m4i  = 


y\  -  j/4 

a;i  —  0:4 


(All) 


J{x-  xkY  +  {y-  ykY  + 


(A12) 


ek  =  {x-  Xkf'  + 


(A.13) 


hk^{x-  Xk){y  -  yk) 


(A.14) 


where  fc  =  1, 2, 3, 4. 


The  velocity  components,  based  on  the  results  from  Hess  and  Smith  [12]  are 


y2-yi  ri+r2-di2  I  yz-y2  r2+r2~d2Z 
d\2  ^  7’i4-r2+di2  '  CZ23  ^  r2+r3+d23 

A-K  I  Vi-yz  In  r3+r4-d34  i  Vi-VA  l„  T^+rx-dji 
dzA  7'3+’'4+«i34  c^41  r4+r  1+^41 


(A.15) 


XX-X2  1^  ri+r2-di2  I  a;2-a^3  1^  T’2+r3-d23 

~  ^12  r\+r2+di2  ^  d23  ^  r2+r3+d23 

47r  I  X3.-a.4  1  r3+r4-ci34  r  xa-Xi  i  rj+n-dn 

d34  ^  r3+r4+<i34  d4i  ^  r*4+r  1+^41 


(AI6) 


Tn23^3—h3 
.  zrs 


^  ^  tan-i  -  tan-i  )  +  tan-^  -  tan-^ 

+tan-i  _  tan-i  +  tan-^  -  tan"! 

(A.17) 

The  u  and  v  components  of  the  velocity  are  defined  everywhere  but  at  the  edges  of  the 
quadrilateral,  where  they  become  infinite.  If  the  influence  of  the  element  on  itself  is  sought, 
then  these  velocity  components  approach  zero  near  the  centroid.  The  jump  in  the  normal 
velocity  component  as  z  — »•  0  inside  the  quadrilateral  is 
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(A.18) 


When  the  point  on  interest  P  {z  =  constant)  lies  outside  of  the  quadrilateral,  then 

u;(>2:  =  0±)  =  0.  (A. 19) 

A  subroutine  was  written  to  incorporate  code  to  handle  nonplanar  panels.  Basically  any  rec¬ 
tilinear  panel  is  input  into  this  subroutine  which  then  adjusts  the  panel  dimensions  to  correctly 
represent  it  as  a  quadrilateral  panel  with  a  x  — y  planar  surface  on  a  constant  z~  surface.  This 
subroutine  also  calculates  the  center  of  the  rectilinear  panel  which  maintains  the  orientation 
between  the  local  (panel)  and  global  (blade  geometry)  coordinate  systems. 

A.2  Fcir  Field  Formulation 

For  improved  computational  efficiency,  when  the  point  of  interest  P  is  far  from  the  center 
of  the  element  (xo^yo,z  =  constant)  then  the  influence  of  the  quadrilateral  element  with  an 
area  of  A  can  be  approximated  by  a  point  source.  When  the  term  “far”  is  invoked,  that  is  a 
distance  more  than  five  average  panel  characteristic  lengths  (longest  distance  between  any  two 
endpoints) ,  then  a  simplified  approximation  is  used.  The  point  source  influence  for  the  velocity 
potential  becomes 


A'k^{x  -  xq)'^  +  {y-  yo)^  + 
The  velocity  components  of  this  source  element  are 

^ _ (taA{x  -  a;o) _ 

47r  [(x  -  xq)^  +  {y  -  yo)^  + 

^  ^ _ o~AA(y  -  yo) _ 

47r  [(a;  -  xo)^  +  (y  -  yo)^  + 

_ (^aA{z  -  ^o) _ 

47r  [(x  -  xo)^  +  (y  -  yo)^  + 


(A.20) 


(A.21) 


(A.22) 


(A.23) 
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An  algorithm  for  calculating  the  influence  of  a  quadrilateral  constant  strength  source  element 
(Program  11)  was  used  [11]. 
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Appendix  B 


Infinite-Pitch  Stator  PBD-I4  Input 
File 
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PBD-14.3  admin  file  for  Infinite-pitch  stator  DSH  6/5/00 


stator. bsn 

File  name  for  blade  b-spline  net 

samp.vel 

File  name  for  wake  field 

1 

15 

15 

nblade,  nkey,  mkey 

1 

1 

ispn  (0==uniform,l=cos),iffixlat 

14 

1,2, 

40 

mctrp,mc(n) 

7 

0.0 

7 

0.0 

ihub,  hgap,  iduct,  dgap 

0 

0 

IHUB  SUBLAY,  IDUC  SUBLAY 

10 

8 

HDWAK,  NTWAKE 

0.84 

300 

Cq,  MXITER 

0.3 

0.3 

OVHANG(l),  OVHANG(3) 

10 

-10 

0 

8 

nx,  ngcoeff,  mltype,  mthick 

4 

IMODE 

0 

NXIMAX 

5 

.001 

.02 

5. 

1 

niter,  tweak,  bulge,  radwgt,  nufix 

1 

.08 

nplot,  hubshk 

4 

2 

NOPT,  NBLK 

999. 

1.0 

1.5 

.05 

ADVCO,  XULT,  XFINAL,  DTPROP 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

:G 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

.95 

1.0 

:r/R 

.015 

.015 

.015 

.015 

.015 

.015 

.015 

.015 

.015 

.015 

:t/D 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

:Cd 

Table  B.l:  Sample  input  file  for  PBD-14. 
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